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Abstract. We introduce and discuss a simple Hamiltonian dynamical system, 
interpretable as a 3-body problem in the (complex) plane and providing the 
prototype of a mechanism explaining the transition from regular to irregular 
motions as travel on Riemann surfaces. The interest of this phenomenology 
— illustrating the onset in a deterministic context of irregular motions - is 
underlined by its generality, suggesting its eventual relevance to understand 
natural phenomena and experimental investigations. Here only some of our main 
findings are reported, without detailing their proofs: a more complete presentation 
will be published elsewhere. 
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Purpose and scope of this paper is to introduce and discuss a simple Hamiltonian 
dynamical system describing the motion of 3 particles in the (complex) plane. This 
3-body problem is the prototype of a class of models [1-3] that feature a transition 
from very simple (even isochronous) to quite complicated motions characterized 
by a sensitive dependence both on the initial data and the parameters ("coupling 
constants") of the model. This transition can be explained as travel on Ricmann 
surfaces. The interest of this phenomenology - illustrating the onset in a deterministic 
context of irregular motions - is underlined by its generality [1,3], suggesting its 
eventual relevance to understand natural phenomena and experimental investigations. 
The novelty of the model treated herein is that it allows a quite explicit mathematical 
treatment. Here only some of our main findings are reported, without detailing their 
proofs: a more complete presentation will be published elsewhere [4] . 

The idea that the integrable or nonintegrable character of a dynamical system is 
closely related to the analytic structure of its solutions as functions of the independent 
variable ( "time" , but considered as a complex variable) goes back to such eminent 
mathematicians as Carl Jacobi, Sophia Kowalewskaya, Henri Poincare, Paul Painleve 
and his school. Some of us heard illuminating discussions of this notion by Martin 
Kruskal, whose main ideas - a synthetic if overly terse rendition of which might be 
the statement that integrability is compatible with the presence of multivaluedness 
but only provided this is not excessive - can be gleaned from some papers written 
by himself and some of his collaborators [5] , or by others who performed theoretical 
and numerical investigations motivated by his ideas [6]. The results presented below 
constitute progress along this line of thinking. For a more detailed analysis we refer 
the interested reader to the papers where more complete versions are presented of our 
findings [4]. 



2. Results 

The model wc introduce and discuss in this paper is characterized by the following 
equations of motion: 

, Sn+2 , 9n+l 

Zn = -lU}Zn-\ 1 . (1) 

Zn -^n+l Zn Zn-\-2 

Notation: here and hereafter indices such as n, m range from 1 to 3 and arc 
defined mod(3); superimposed dots indicate differentiations with respect to the real 
independent time variable t; the dependent variables Zn = Zn{t) are complex, and 
indicate the positions of 3 point "particles" moving in the complex z-planc: i = i/ — 1 
is the imaginary unit; the parameter lv is positive, and it sets the time scale via the 
basic period 

T^l-. (2) 

the 3 quantities gn are arbitrary coupling constants, but in this paper we restrict 
consideration to the case in which they are all real and moreover satisfy the 
"semisymmetrical" restriction 

gi= 92= g , 93 = f , (3) 
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entailing that the two particles with labels 1 and 2 are equal, while particle 3 
is different. More special cases are the "fully symmetrical", or"integrable" , one 
characterized by the equality of all 3 coupling constants, 

/ = 3 , 51 = .92 = 53 = .9 , (4) 
and the "two-body" one, with only one nonvanishing coupling constant, 

51 = 52 = 5 = 0, 53 = / ^ . (5a) 
In this latter case clearly 

23(i) = ;23(0) exp(-icji) (56) 
and the remaining two-body problem is easily solvable, 



I [Zl(0) + Z2(0)] 



s{t) = exp(— icjt) 

exp(2itjt) - 1-1 1/2 



-(-r{^[^i(o)-.2(o)p + /^^^^^^^ 



1,2. (5c) 



The justification for labelling the fully symmetrical case l@J as "integrable" will be 
clear from the following (or see Section 2.3.4.1 of [3]). The treatment of the more 
general case with 3 different coupling constants is outlined in [4]. 

Note that the equations of motion are of "Archimedian" , rather than 
"Newtonian", type, inasmuch as they imply that the "velocities", rather than the 
"accelerations", are determined by the "forces". These equations of motion are 
Hamiltonian, indeed they follow in the standard manner from the Hamiltonian 
function 



n=l 



Pn+1 - Pn+2 
-1 W Z„ p„ + 5„ 



(6) 



And they can be reformulated [4] as, still Hamiltonian, real (and covariant, even 
rotation-invariant) equations describing the motion of three point particles in the 
(real) horizontal plane. 

The following qualitative analysis (confirmed by our quantitative findings, see 
below) is useful to get a first idea of the nature of the motions entailed by our model. 
For large values of (the modulus of) 2„ the "two-body forces" represented by the last 
two terms in the right-hand side of become negligible with respect to the "one- 
body (linear) force" represented by the first term, hence in this regime i;„ ~ — i w z„ 
entailing z„(t) « const exp{—iujt). One thereby infers that, when a particle strays far 
away from the origin in the complex z-plane, it tends to rotate (clockwise, with period 
2 T) on a circle: hence the first qualitative conclusion that all motions are confined. 
Secondly, the two-body forces cause a singularity whenever there is a collision of two 
(or all three) of the particles, and become dominant whenever two particles get very 
close to each other, namely in the case of near misses. But if the three particles move 
aperiodically in a confined region (near the origin) of the complex z-plane, an infinity 
of near misses shall indeed occur. And since the outcome of a near miss is generally 
quite different (whenever the two particles involved in it are different) depending 
on which side the particles slide past each other - and this, especially in the case 
of very close near misses, depends sensitively on the initial data of the trajectories 
under consideration - we see here a mechanism causing a sensitive dependence of the 
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time evolution on its initial data. This suggests that our model (Q, in spite of its 
simplicity, might also support quite complicated motions, possibly even displaying 
an "unpredictable" evolution in spite of its deterministic character. This hunch is 
confirmed by the results reported below. 

To investigate the dynamics of our "physical" model it is convenient to 
introduce an "auxiliary" model, obtained from it via the following change of dependent 
and independent variables: 

Zn it) = exp {-I LO t) Cn (t) , T (t) = ^. ' . (7) 

Note that initially the coordinates and c^^ coincide: 

2„(0) = Cn(0) . (8) 

The equations of motion of the auxiliary model follow immediately from ^ via 
(or, even more directly, by noting that, for a; = 0, r = i and Zn{t) = Cn (■'')): 

>/ 9n+2 , gn+l /„x 

C„ - 7 7 1- 7 7 — ■ (^) 

Here of course the appended prime denotes differentiation with respect to the 
{complex) variable r. 

The definition of r (t) implies that as the (real) time variable t evolves onwards 
from t = 0, the complex variable r travels round and round, making a full tour 
(counterclockwise) in every time interval T, on the circle C the diameter of which, of 
length d = 1 / Lo, lies on the imaginary axis in the complex r-plane, with one end at 
the origin, t = 0, and the other at r = i / w (draw this circle!). Hence these relations, 
l|7|l. entail that if C„ (r) is holomorphic as a function of the complex variable t in the 
closed disk D encircled by the circle C, the corresponding function z„ (t) is periodic 
in the real variable t with period 2 T (indeed antiperiodic with period T) : 

z„(t + T) = -2„(i) , z„(t + 2r) = z„(i) . (10) 

But it is easy to prove [4] that the solution C„ (r) of ^ is holomorphic (at least) in the 
circular disk Dq centered at the origin of the complex r-plane and having the radius 

2 



min |C„(0)-C™(0)| 

n,m— 1,2, J; m^n 

128 max 

n— 1,2,3 



(11) 



One may therefore conclude that our physical system |^ is isochronous with period 
2 T, see (jSJ . Indeed an isochronous system is characterized by the property to possess 
one or more open sectors of its phase space, each having of course full dimensionality, 
such that all motions in each of them are completely periodic with the same fixed 
period (the periods may be different in these different sectors of phase space, but must 
be fixed, i e. independent of the initial data, within each of these sectors): and in our 
case clearly (at least) all the motions characterized by initial data z„ (0) such that 



/ max \gn\ 

\z^{Q)-z^m>^^\r^^ — (12) 



n,m— 1,2,3; m^n V ^ 

are completely periodic with period 2T, see Hl()(l . since this inequality, implying (via 
(jHJ and Hll|l ) r > d, entails that Cn{T) is holomorphic (at least) in a disk Dq that 
includes, in the complex r-plane, the disk D. 
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This argument is a first demonstration of the usefulness of the "trick" Q, 
associating the auxihary system ^ to our physical system More generally, 

this relationship |(7J) allows to infer the main characteristics of the time evolution 
of the solutions Zn (t) of our physical system from the analyticity properties of 
the corresponding solutions Cn (t) of the auxiliary system indeed the evolution of 
Zn{t) as the time t increases from the initial value t = is generally related via 
to the values taken by Cn (t) when r rotates (counterclockwise, with period T) on the 
circle C in the complex r-plane and correspondingly Cn{T) travels on the Riemann 
surface associated to its analytic structure as a function of the complex variable r. 
Suppose for instance that the only singularities of Cn{T) in the finite part of the 
complex T-plane are square-root branch points, as it is indeed the case for our model 
Q at least for a range of values of the ratio of the coupling constants / and g, see 
lpj)l and below.l Then the isochronous regime corresponds to initial data such that 
the corresponding solution C„ (r) has no branch points inside the circle C on the main 
sheet of its Riemann surface (i. e. that characterized by the initial data). Moreover, 
if there is a finite (nonvanishing) number of branch points inside the circle C on the 
main sheet of the Riemann surface of Cn (t) , and a finite number of branch points 
inside the circle C on all the sheets that are accessed by traveling on the Riemann 
surface round and round on the circle C, then clearly the corresponding solution Zn {t) 
is still a completely periodic function of the time t, but now its period is a finite integer 
multiple jT of the basic period T, the value of j depending of course on the number of 
sheets that get visited along this travel before returning to the main sheet. Hence, in 
particular, whenever the total number q of (square-root) branch points of the solution 
Cn{T) of the auxiliary problem Q is finite, the corresponding solution z„(t) of our 
physical model is completely periodic, although possibly with a very large period 
(if q is very large) the value of which may depend, possibly quite sensitively, on the 
initial data. On the other hand if the number of {square-root) branch points possessed 
by the generic solution CnC^) is infinite, and the Riemann surface associated with 
the function CniT) has an infinite number of sheets (as it can happen in our case, 
see below), then it is possible that, as r goes round and round on the circle C, the 
corresponding value of ^„ (t) travels on this Riemann surface without ever returning 
to its main sheet, entailing that the time evolution of the corresponding function z„(i) 
is aperiodic, and that it depends sensitively on the initial data inasmuch as these data 
characterize the positions of the branch points hence the structure of the Riemann 
surface. 

This terse analysis entails an important distinction among all these (square-root) 
branch points: the "active" branch-points are those located inside the circle C on 
sheets of the Riemann surface accessed - when starting from the main sheet - by 
traveling round and round on that circle, so that they do affect the subsequent sequence 
of sheets that get visited; while the "inactive" branch points are, of course, those that 
fall outside the circle C, as well as those that are located inside the circle C but on 
sheets of the Riemann surface that do not get visited while traveling round and round 
on that circle (starting from the main sheet) and that therefore do not influence 
the time-evolution of the corresponding solution of our physical system (^. This 

X The nature of these singularities is generally independent from the particular solution under 
consideration and can be easily ascertained [4] via local analyses d la Painleve of the generic solution 
of the equations of motion while the number and especially the locations of these singularities 
depend on the specific solution under consideration and their identification requires a more detailed 
knowledge than can be obtained by a local analysis d la Painleve. 
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distinction is of course influenced by the initial data of the problem, that characterize 
the initial pattern of branch points; clearly it is not just a "local" characteristic of 
each branch point depending only on its position (for instance, inside or outside the 
circle C): it depends on the overall structure of the Rieniann surface, for instance if 
there is no branch point on its main sheet - that containing the point of departure of 
the travel round and round on the circle C - then clearly all the other branch points 
are inactive, irrespective of their location. 

Let us also emphasize that, whenever an active branch point is quite close to the 
circle C, it corresponds to a near miss involving two particles of our physical model 

at which these two particles scatter against each other almost at right angles 
(corresponding to the square-root nature of the branch point). The difference between 
the cases in which such a branch point falls just inside respectively just outside the 
circle C corresponds to a near miss in which the two particles slide past each other 
on one side respectively on the other (see figure ^ , and this makes a substantial 
difference as regards the subsequent evolution of our 3-body system (unless the two 
particles are equal). The closer the near miss, the more significant this effect is, and 




t - plane plane 



Figure 1. Scattering of two bodies in the three-body problem Q corresponding 
to a near miss. The two outcomes originate from two sets of initial data close to 
each other such that a square-root branch point rj, falls on different sides of the 
circle C. 

the more sensitive it is on the initial data, a tiny change of which can move the 
relevant branch point from one side to the other of the circumference of the circle 
C and correspondingly drastically affect the outcome of the near miss. This is the 
mechanism that accounts for the fact that, when the initial data are in certain sectors of 
their phase space (of course quite different from that characterized by the inequalities 
H12|l '). the resulting motion of the physical 3-body problem ^ is aperiodic, indeed 
nontrivially so: in such cases (as we show below) the aperiodicity is indeed associated 
with the coming into play of an infinite number of (square-root) branch points of 
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the corresponding solution of the auxiliary problem ^ and correspondingly with an 
infinite number of near misses experienced by the particles throughout their time 
evolution, this phenomenology being clearly characterized by a sensitive dependence 
on the initial data. 

This mechanism to explain the transition from regular to irregular motions - 
and in particular from an isochronous regime to one featuring unpredictable aspects 
- was already discussed [2] in the context of certain many-body models somewhat 
analogous to that studied herein. But those treatments were limited to providing a 
qualitative analysis such as that presented above and to ascertaining its congruence 
with numerical solutions of these models. The interest of the simpler model introduced 
and discussed herein is to allow a detailed, quantitative understanding of this 
phenomenology. This is based on the following explicit solution of our model QJ, 
obtained [4] via the auxiliary problem © : 



Zs{t) = Z exp {-iut) - - 



1 ff + SgV^' _^ , „. ,,,1/2 

2 ( g-^ j [1 + ?7 exp(-2iwt)] ' • 

■{[w{t)]'/^~{~r [U^^~3w{t)]'^'} , s = l,2, (13a) 

Z3{t) = Z exp{-iujt)+ (^-^] ^ [1 + r; exp(-2ic^t)]'/' [z«(t)]'/' . (13&) 

\ biuj J 

Here the function is defined via the relation 

wit)^wm], (14) 

with 

C{t) ^ R [ri + e^p{2iujt)] = ^ + R e^p{2iujt) , (15) 

and w{^) implicitly defined by the nondifferential equation 

(tu-l)^-^ u;-'' =^ . (16) 

The parameter fj, is defined in terms of the coupling constants g and /, see ©, as 
follows: 

and in fl3V (|15|) the three constants Z, R, and rj (or ^) are defined in terms of the 3 
initial data z„(0) as follows: 

^ Zl(0)+.22(0) + Z3(0) 



R 



3 

3(/ + 85) ' ' ^"-^ 



2iuj [2z3(0)-zi(0)-Z2(0)]^ 



(18a) 
(18&) 



w(0) 

e-i?77, (18c) 

i { [Zl (0) - Z2 (0)] ' + [Z2 (0) - Z3 (0)] V [Z3 (0) - Zl (0)] ' } 



3 {f + 2g) 



(18d) 



I 2 



2m[2z3(0)-^i(0)-^2(0)]- . 

wlO = ^ ^ ^ . (18e) 

[Zl (0) - Z2 (0)] V [Z2 (0) - Z3 (O)] V [Z3 (O) - zi (0)] ' 
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Note that the constant Z is the initial value of the center of mass of the system, 
and indeed the first term in the right-hand side of the solution (13) represents the 
motion of the center of mass of the system: just a circular motion around the origin, 
with a constant velocity entailing a period 2 T. Since the rest of the motion is 
independent of the behavior of the center of mass, in the study of this model attention 
can be restricted without significant loss of generality to the case when the center of 
mass does not move, Z = 0. 

The nontrivial aspects of the motion are encoded in the time evolution of the 
function 'w{t), see (13) and (|14|) : let us emphasize in this connection that the dependent 
variable w{^) is that solution of the nondifferential equation (|16|l uniquely identified 
by continuity, as the time t unfolds, hence as the variable ^ = ^ (t) goes round and 
round, in the complex ^-plane, on the circle S with center ^ and radius \R\ (see H15|l 'l. 
from the initial datum assigned at i = 0, 



see (|18e|) . This specification of the initial value ■!i7(0) is relevant, because generally 
the nondifferential equation has more than a single solution, in fact possibly an 
infinity of solutions, see below. 

It is clear from (13) that the time evolution of the solution z„(t) of our model Q 
is mainly determined by the time evolution of the function w{t). Indeed, 

1/2 

(i) The factor [?/ exp (—21 wt) — 1] displays a quite simple time evolution, periodic 
with period T if jr^l < 1 and antiperiodic with period T hence periodic with period 
2r if |7/| > 1. 

1/2 

(ii) If 'w{t) is periodic with period T, its square root [w{t)] , appearing in the right- 
hand side of the solution formulas (13), is clearly as well periodic with period T 
or antiperiodic with period T hence periodic with period 2 T depending whether 
the closed trajectory of ■w{t) in the complex w-plane does not or does enclose the 
(branch) point w = 0. 

(iii) Likewise the square root [I2fi — 5w (t)] (see (I13afl ^ is also periodic with period 
T or antiperiodic with period T hence periodic with period 2 T depending whether 
the closed trajectory of w{t) in the complex w-plane does not or does enclose the 
(branch) point w = Afx (but note that a change of sign of this square root only 
entails an exchange between the two equal particles 1 and 2). 

(iv) In conclusion one sees that - provided one considers particles 1 and 2 as 
indistinguishable - then, if the time evolution of 'w{t) is periodic with period 
T, w(i + T) — w(t), the physical motion of the 3 particles Zn(t) is also completely 
periodic either with the same period T or with period 2 T, provided T is an integer 
multiple of T; 

(v) Finally, if the motion of w^t) is not periodic then clearly the functions Zn{t) are 
also not periodic. 

Hereafter we only discuss the time evolution of the function wit)] actually, as 
explained below, in this paper we limit our consideration to discussing the motion of 
a generic solution 'w{t) — w [^{t)] of the nondifferential equation (|16() . Moreover, we 
consider only generic solutions of Q , namely those characterized by initial data that 
exclude one of the following special outcomes: 

(a) w (t) takes, at some (real) time ta, the value w{ta) = 4:fj, entailing a pair collision 
of the 2 equal particles occurring at this time, Zi{ta) = Z2(ta)- 




(19) 
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(b) w (t) takes, at some (real) time tf,, the value w(t{,) = /i entailing a pair collision 
of the different particle with one of the 2 equal particles occurring at this time, 
ziih) = z?,{h) or Z2(4) = z^ih)- 

(c) The constant 77 has unit modulus, It^I = 1, i. e. rj — exp(2icjtc) with tc real 
(and of course defined mod T) which entails a triple collision of the 3 particles 
occurring at the time tc, zi{tc) = Z2{tc) = zj,{tc) . 

(c') w (t) vanishes at some (real) time t^ w{tc) = (but, as our notation suggests, 
this case (c') is just a subcase of (c), although this is not immediately obvious 
from (|13all but requires using also IjlSfl and (|16|l ). 

The initial data that give rise to solutions having one of these singularities form a set 
of null measure. It can be easily seen that these singular solutions Zn{t) of our physical 
problem correspond via Q to special solutions Cn (t) of our auxiliary problem (|5J 
possessing a branch point that sits exactly on the circle C in the complex r-plane: 
more precisely, 

(a) a sguare-root branch point featured by Ci(''') ^^^id C2(''') but not by Csl''"): 

(b) a square-root branch point featured by all 3 functions Cti(''')i 

(c) a branch point featured by all 3 functions Cn(''") the nature of which depends on 
the parameter ^. 

As mentioned above, in this paper we confine our treatment to discussing the time 
evolution of a generic root w{t) of the nondifferential equation H16|l with H15() . and 
in particular to identifying for which initial data its time evolution is periodic, and in 
such a case what the period is. Remarkably we find out that, for (arbitrarily) given 
initial data, all these roots have at most three different periods (one of which might 
be infinite, signifying an aperiodic motion); periods which we are able to determine 
explicitly (although the relevant formulas have some nontrivial, even "chaotic" , 
aspects, in a sense that is made explicit below). The question of identifying, among 
all the roots 'Wj{t) of this nondifferential equation H16|l. the "physical" one 'w{t) i. e. 
the one that evolves from the initial datum H18ell , and in particular of specifying the 
character of its time evolution among the (at most 3) alternatives discussed below, 
is a technically demanding job the solution of which shall be reported in [4]. Let us 
re-emphasize that the time evolution of w{t) = w [£,{t)] coincides with the evolution of 
a generic root w (^) of H16|) as the independent variable ^ travels (making a complete 
counterclockwise tour in the complex ^-plane in every time interval T) on the circle 
S with center ^ and radius \R\, see ()15f) . and correspondingly the dependent variable 
w{£,) travels on its Riemann surface. Note that this Riemann surface is completely 
defined by the single parameter /i, see (|17|l and (I16|l . while the circle S is defined by 
the initial data of the problem, see l|15|l with (18). 

What therefore remains to be discussed is the analytic structure of the multivalued 
function w (^) defined implicitly by the nondifferential equation or, equivalently 
but more directly, the time dependence of the corresponding function w (t) = w (t)]. 
To begin with, we consider the case in which the parameter /i is rational, 

M - - , (20) 
9 

with p and q coprime integers and q positive, q > 0. The extension of the results to the 
case of irrational fj, is made subsequently; although, to avoid repetitions, we present 
below some results in a manner already appropriate to include also the more general 
case with n real. 
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In the rational case (|20|l the nondijferential equation that determines the 
"dependent variable" w (t) in terms of the "independent variable" t becomes 
polynomial, and takes one of the following 3 forms depending on the value of the 
parameter yit, see ^U^: 

{w-l)P~'^ ^[^+Rexp{2iujt)]'^ wP, if^>l, (21a) 

[^+Rexp{2iujt)]'' (w-l)"^" w" ifO<Ai<l, (21&) 

[C + i?exp(2iwt)]'' (i&-l)''+IPl =i&IPl, ifAi<0. (21c) 
The above expressions are polynomials (in the dependent variable w) of degree J: 

{p, if > 1; 

q, if < < 1; (22) 

q+\p\, if /^<0. 

As for the boundaries of these 3 cases, let us recall that ^ = 1 corresponds, via H17|l . 
to g = 0, namely, see (O, to the trivially solvable two-body case, see (5), while ^ — 
respectively fj, = oo correspond, via (|17|l . to / + 2 g = respectively to / + 8 g = and 
require a separate treatment, for which the interested reader is referred to [4]. Clearly 
the third case (/i < 0) becomes identical to the first (/i > 1) via the replacement 

wi-^l — w, S.'-^—S.i i? 1-^ — i?, —p^P'-q, q ^ p i-^ p (23) 
without modifying q; therefore in the following, without loss of generality, we often 
forsake a separate discussion of this third case. 

Clearly the factor + i? exp(2iwt)] , which carries all the time dependence 
in these polynomial equations, is periodic in t with period T, see (except for the 
special initial conditions entailing ^ = 0, in which case this factor is instead periodic 
with the shorter period T / q) (for simplicity we continue to pursue our policy to 
consider only the generic case when this does not happen, referring the interested 
reader to [4] for a more complete treatment). At issue is the behaviour of the J roots 
Wj(t) of our polynomial equation (21) whose coefficients evolve in time periodically 
with period T . Let us indicate with W{t) = {'Wj{t)\ j — 1, J} the (unordered) set 
of these J roots. Obviously W{t) is periodic with period T, W{t + T) = W{t): after 
one period T the polynomial equation is unchanged, hence the set of its J roots is 
as well unchanged. But that does not imply that if one follows the time evolution 
of these J roots, each of them will return to its own initial value after one period, 
Wj(T) = Wj{0), j = 1, J. This outcome will indeed obtain for the open domain of 
initial data of our problem that corresponds to the basic isochronous regime, see (|10|) : 
but it does not happen for other initial data, in which cases for instance a generic 
root, say wj-^ (t) = w{t), may after one period land in the initial position of a different 
root, say w{T) = ^^2(0), and then after one more period end up in the initial position 
of yet another root, w(2T) = Wj^ijd), and so on. Eventually, of course, after a time 
T = jT which is a finite integer multiple j of the basic period T, with 1 < j < J, 
the generic root w{t) shall necessarily return to its initial position, w (^T^ = w (0) , 
entailing that its evolution as a function of the time t is periodic with this period T, 
w(t + f^=w{t + jT)=w{t) . (24) 

This discussion clearly implies (via (|13Qfl ^ that, in the case now under consideration 
(with a rational value of the parameter /i, see H2()|l ^. all solutions of our physical 
problem with (jSjl are completely periodic with a period which is either 

f = jT with 1 < j < J (25) 
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and J defined by (|22|) . or it is 2f (see the discussion above following (|19|l 'l. The 
remaining, crucial question is: how does the value of the integer j (which might be 
quite large if J is quite large) depend on the initial data of our problem? In this paper 
we outline how to calculate, for given initial data, all the possible periods of the J 
roots Wj{t) of (21), and we display formulas providing (at most) 3 alternative values 
for these periods; as already mentioned above, the explanation of how to identify 
which one of these 3 periods corresponds to that of the "physical" root w{t) entails a 
more detailed treatment which is reported in [4]. 

But before doing so let us emphasize that via this discussion the time evolution of 
our original 3-body problem - describing the time evolution of the three points Zn(t) 
in the complex z-plane - has been related to the time evolution of the J roots Wj (t) of 
(21) in the complex u;-plane, and in particular to the way they get permuted among 
themselves over the time evolution after each period T. As we will explain below, the 
possible complications of the motions of our physical model ^ are thereby related to 
the mechanisms at play to permute these roots among themselves when one watches 
their time evolution at subsequent intervals T, 2 T, 3 T and so on. And, in this context, 
it is significant to note that, whenever fi is irrational, one is in fact dealing with the 
dynamics of an infinite number of roots. This suggests that whenever the number 
J of roots is large, and even more so when /i is irrational (entailing J = oo), the 
time evolution of our physical model ^ might be quite complicated, perhaps calling 
into play the theoretical tools of statistical mechanics rather than those of few-body 
dynamics (but we postpone such excursions to future publications). 

As entailed by our discussion above, the issue of determining the value 
of the integer j is tantamount to understanding the structure of the Riemann 
surface associated with the function w(^) of the complex variable ^ defined by the 
nondifferential equation H16|l . whose different sheets correspond of course to the 
different roots of this nondifferential equation. In the rational case ()20|l this equation 
is in fact polynomial of degree J in the dependent variable w. Specifically, what must 
be ascertained is the number of sheets of this Riemann surface that are accessed by 
w;(^) when the independent variable ^ travels in the complex ^-plane round and round 
on the circle S, whose center ^ and radius \R\ depend on the initial data of our physical 
problem, see (I15|) and (18). To this end one must gain and use a detailed understanding 
of the structure of this Riemann surface. Again, we refer for the details of this analysis 
to [4] . Suffice here to report the following basic information on the analytic structure 
of the function w (^) , referring to the general case with real fi {rational or irrational). 

The nondifferential equation (|16|) defines a J-sheeted covering of the complex ^- 
plane of genus zero (of course J = cx) is ^ is irrational). The function w{^) defined 
implicitly by this equation features square-root branch points located on a circle B 
centered at the origin of the complex ^-plane: 

6=^f^ -?'6exp(2^i^fc) , A: = 1,2,3,... , (26a) 



^b = Cb =n cxp 



rb = {^,-l)-' . (26c) 

In the last, (|26c|) . of these formulas it is understood for definiteness that the principal 
determination is taken of the fi-th power appearing in the right-hand side. The first of 
these formulas, H26a|) . shows clearly that the number of these branch points is infinite 



, 2 np k 
1 



fc = l,2,...,g, (26&) 



fi-l 
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if the parameter fi is irrational, and that they then sit densely on the circle B in the 
complex ^-plane centered at the origin and having radius |rf,|, see (|26ef) . Note that 
this entails that the generic point on the circle B is not a branch point (just as a 
generic real number is not rational); but every generic point on the circle B has, if fj, 
is irrational, some branch point (in fact, an infinity of branch points!) arbitrarily close 
to it (just as every generic real number has an infinity of rational numbers arbitrarily 
close to it). It is also important to realize that these branch points are generally on 
different sheets of the Riemann surface associated with the function w {£_): hence, they 
are dense if one considers the circle B in the complex ^-plane, but they are not dense 
if one considers these branch points on the Riemann surface itself. As for the second 
of these formulas, Ij266|) . it is instead appropriate to the case in which the parameter 
fx is rational, see ()20|l , in which case the branch points sit again on the circle B in the 
complex ^-plane, but there are only a finite number, q, of them (and note that the 
factor p appearing in the argument of the exponential in the right-hand side of this 
formula, (j26&|) . is only relevant to characterize how the branch points ^^''^ are labeled 
via the index k). Both in the irrational and in the rational case at all these branch 
points the nondifferential equation has a double root which takes the same value 
ii;(^f,) = fi. Note that this entails that circling around such a branch point in the 
complex ^-plane corresponds to permuting 2 of the J roots Wj (^) among themselves. 

In addition, the function w{^) possesses branch points at ^ = and at ^ = oo, the 
order of which depends on the value of ft, and is rational if the number fj, is rational 
(see below). 

The branch point at ^ = cx3 has, if /i > 1, exponent — ^ {— — ^ in the rational 
case), 

^«oo, w{£_) ^ aC^/^" ^0 , a^" = -exp{-iTrfi) , fi > , (27) 
while it has instead two different exponents if < /i < 1: 

(i) the exponent — ^ ^^^'^ rational case) as given by the preceding formula 
(entaihng w « 0) 

(ii) the exponent ^jr^ (= in the rational case) as given by the following formula 
(entailing w « 1), 

C«oo, w(0 « l + ar^/'^""' « 1 , 0^-^ = 1, 0<fL<l. (28) 
The branch point at ^ = is, if /i > 1, of exponent (= in the rational 
case), 

^«0, w{0~l + a^,^/^^'-^'> ^1 , a^-i = l, fi > 1 , (29) 

(note the formal analogy of this formula with the previous one, (jSHl), and it is instead 
absent if < /i < 1, so that in this second case the only branch points in the finite 
part of the complex ^-plane are those of square-root type, see (|26a|l . This is the 
main cause of the difference between the results, see below, for this case (0 < ^ < 1) 
from those for the other two cases (/i > 1 and /i < 0), which are on the other hand 
essentially equivalent among each other being related by the transformation /i i— s- I — fi, 
w{£,) I— > 1 — w {—£,) , see (|16|) (so that without loss of generality we often forsake an 
explicit discussion of the case with fi < 0). 

Note that these results entail that, in the rational case, see (|20|) . making a circle 
around the branch point at ^ = 0, in the p > q {i. e. /i > 1) case when this 
branch point is present, causes a cyclic permutation of p — g of the p roots wj: this 
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is particularly evident if one imagines to travel full circle around the branch point 
at ^ = in its immediate vicinity, since for ^ « the p roots Wj of (|16|l are clearly 
divided into two sets, a first set oi p — q roots, disposed equispaced on a circle of small 
radius (« ICI'^'^''^ centered at w = 1 in the complex w-plane, which then undergo a 
cyclic permutation among themselves, and a second set of q roots, disposed equispaced 
on a circle in the complex w-plane centered at the origin and having a large radius 
(« each of which after the operation returns instead to its original position. 

The permutation experienced by the J roots lij (t) due to a sequence of pairwise 
exchanges of roots - which take place whenever square-root branch points are included 
inside the circle S traveled by the point ^{t) - causes a reshuffling of the roots which 
is nontrivial inasmuch as it depends on how many and on which pairs of roots get 
sequentially exchanged over each period, as determined by the number and identity 
of square-root branch points enclosed inside the circle S traveled by ^(t) and by 
the detailed structure of the Riemann surface associated with these branch points, 
in particular on which sheets of this Riemann surface the relevant branch points 
are located. The reshuffling encompasses more roots when a second mechanism is 
simultaneously at play, i. e. that producing a cyclic permutation oi p — q roots 
(especially, of course, when p — q is large) over each period, as caused by the presence 
of the branch point at ^ = 0: this second mechanism exists only if > 1 (or /i < 0, 
entailing the exchange —p ^ p — q), and provided the circle S traveled by ^(t) does 
include the point ^ = 0. This phenomenology causes the possible periods T — jT 
of the time evolution of the generic root w{t) to depend on the initial data, but 
remarkably we will see that, for given initial data, there are (at most) only 3 possible 
values of these periods, and that they can be given explicitly in terms of the initial 
data and of the two numbers p and q, see (|2U|I . Indeed we show below that analogous 
results can as well be given in the case when /i is irrational, in spite of the fact that 
the dependence on the initial data may then be quite sensitive. 

The correspondence of the analysis, given here in terms of the dynamics of the 
roots Wj{t), with that in terms of travel on the Riemann surface made above (see 
paragraphs after H12|l ). including the distinction made there about active and inactive 
branch points, should be noted: the active branch points are those that cause a 
reshuffling of roots that involves the "physical" root 'w{t), the inactive ones are those 
that do not cause a reshuffling of roots that involves the "physical" root w(t), either 
because they cause no reshuffling at all being located outside the relevant circle (S in 
the complex ^-plane in the context of the present analysis, C in the complex r-plane 
in the context of the discussion made above when the distinction among active and 
inactive branch points was first introduced), or because they cause a reshuffling which 
however does not involve the physical root w{t). 

We now report our findings concerning the time evolution of a generic root w{t), 
referring at first mainly to the rational case but including immediately results for the 
irrational case whenever it is convenient to do so in order to shorten our presentation. 
The remaining information on the irrational case is provided below (see Proposition 
5). ^ 

First of all it is useful to visualize the two circles B and ^ in the complex ^-plane 
(draw them!): recall that the circle B on which the branch points sit is centered at the 
origin and its radius |rf,| only depends on the parameter ^, see (126 c|l . while both the 
center ^ and the radius \R\ of the circle 5 traveled upon by (,{t), see (|15|l . do depend 
on the initial data, see (18). 



Transition from regular to irregular motions 



14 



Proposition 1. If the circle S is inside the circle B (i. e. |^| + \R\ < \rh\ , see 
(18) and H2t)cll ). and (a) /i is inside the interval < /i < 1 or (b) fi is outside this 
interval (/x > 1 or ^ < 0) and moreover the circle S does not include the origin C = 
(i. e. \R\ < or equivalently |?7| > 1), then j = 1, i. e. the generic solution w{t) is 
periodic with period T, w{t + T) = w{t). This outcome applies equally if ^ is rational 
or irrational. □ 

Remark 1. In all the cases identified in this Proposition 1 there are no branch 
points at all inside the circle S: indeed the outcome detailed by this Proposition 1 
applies in all the cases in which this happens (see our discussion above), even (in 
the case with rational fi) if the circles B and S do cross each other marginally; and 
of course in all these cases all the roots Wj{t) are periodic with period T, and in 
the context of the physical problem the solution is characterized by the simple 
periodicity rule (|10|l . The restriction (|12|l on the initial data is sufficient (but of 
course not necessary) to guarantee that we are in this regime. □ 

Proposition 2. If the circle S is inside the circle B (i. e. |f| + \R\ < \ri,\ ), the 
circle 5 does include the origin ^ = (i. e. \R\ > |^| or equivalently \ri\ < 1), and /x 
is outside the interval < /i < 1 then in the rational case, see (|20|l . 

(a) j = loTj=p— qiffi>l 

(b) j = 1 or j = \p\ if < 0. 

In the irrational case the time evolution of the generic root w^t) is either periodic with 
the basic period T, or quasiperiodic, involving in particular a (nonlinear) superposition 
of two periodic evolutions with two noncongruent periods, specifically (a) with period 
T and ^ if A* > 1, (b) with period T and ^ if ^ < O.D 

Remark 2. In the case identified in this Proposition 2 the only branch point 
inside S is that at ^ = 0, which is indeed only present if /i > 1 or /i < 0. Hence in the 
rational case with p > q {i. e. /i > 1), p — q roots Wj{t) get cyclically exchanged among 
themselves, entailing that the time evolution of each of them has period (p — q) T, 
while the remaining q roots have period T; with an analogous phenomenology in the 
/X < 0. Likewise, when /i is irrational, the periodicity of the time evolution of the 
generic root 'w{t) has period T if the branch point at ^ = is inactive (i. e., it does 
not appear on the sheet on which 15(0) lives), otherwise its time evolution can be 
inferred by replacing ^ with ^(t), see (|15|l . in the formula characterizing the branch 
point at ^ = 0, see (|29|l . Note however that for the special initial data such that the 
two circles B and S are concentric (i. e., 77 = 0) the quasiperiodic time evolution of 
'w{t) is instead periodic (a) with period -^^j if /x > 1, (b) with period if /i < O.D 

Proposition 3. If the circle S is outside the circle B (i. e. \R\ > |^| + \rt,\ ) then 
in the rational case, see (|20|l . 

(a) j j3 if /i > 1, 

(b) j = q-p = q+\p\it fi<Q 

(c) j— porj = q — piiO<ii<l. 

In the irrational case the time evolution of the generic root 'w{t) is quasiperiodic, 
involving a (nonlinear) superposition of two periodic evolutions with two noncongruent 
periods, specifically 

(a) with periods T and ^ if /x > 1, 

(b) with periods T and -j^— if /x < 0, 
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(c) with periods T and ^ or T and if < ^ < l.D 

Remark 3. In aU the cases encompassed in this Proposition 3 the branch points of 
w in the finite part of the complex ^-plane are all inside the circle 5, hence the 
dynamics of the roots Wj(t) can be understood in terms of the branch point of w {^) 
at ^ = oo. Therefore in the rational case with p > q {i. e. /i > 1) all the p roots get 
cyclically exchanged, so that each of them gets back to its original value after a period 
pT; likewise if p < (i. e. /z < 0) all the q + \p\ roots get cyclically exchanged, so 
that each of them gets back to its original value after a period {q + \p\) T. In the other 
rational case, < p < q {i. e. < /i < 1), p roots gets cyclically exchanged among 
themselves, and the remaining q — p roots get cyclically exchanged among themselves, 
so that the generic root w{t) has period pT if it belongs to the first set, and (q — p) T 
if it belongs to the second. And the outcome in the irrational case can as well be 
understood in terms of the exponent of the branch point at ^ = oo, see (|27|I - H28|I .D 

The situation is less straightforward (hence more interesting) if the two circles B 
and S do intersect each other (i. e. |f| — \R\ < \ri,\ < |^| + \R\). Then the parameter 
that plays a crucial role is the number b of square-root branch points, sitting on the 
circle B, that fall inside the circle S. This number b is finite, 1 < b < q (note 
that the case 6 = is taken care of by Proposition 1) only in the rational case, to 
which we restrict consideration in the following Proposition 4 (and we exclude from 
consideration the nongeneric case in which the circle S hits one of the branch points 
sitting on the circle B). Then two or three (but no more!) different alternatives are 
possible for the value of the positive integer j, as detailed below. 

Hereafter the notation \x\ denotes the floor of the real number x, namely the 
largest integer number not larger than x (hence for instance [— 0.3J = — 1, [OJ =0). 

Proposition 4. (i). lip > q {i. e. /i > 1, see (|20|l ') and the origin ^ = is outside 
the circle S (i. e. I77I > 1), then j can take one of the following 3 values: 



i = 1 or i 



1 



p-q 



1 or i 



1 



P-q. 



(30) 



(ii). If p > g (i. e. ^ > 1, see ((201)) and the origin ^ = is inside the circle ^ (i. 
I77I < 1), then j can take one of the following 2 values: 

j = l or j = b + p-q. (31) 



(i'). If p < (i. e. /i < 0, see (|2U|l ') and the origin ^ = is outside the circle ^ (i. 
|?7| > 1), then j can take one of the following 3 values: 



j = 1 or j = 



6-1 
1^ 



1 or j = 



b-1 



(32) 



(ii'). If p < (i. e. /i < 0, see (|20|) ') and the origin ^ = is inside the circle S (i. 
e. I77I < 1), then j can take one of the following 2 values: 

j = l or 3 = b+\p\ . (33) 

(iii). The situation is more intriguing if < p < g (i. e., < < 1, see ((201)). 
Then one must introduce the simple continued fraction expansion of the number 

q 1 1 



Transition from regular to irregular motions 



16 



The fcth convergent Ck of this continued fraction expansion H34|l . and in particular its 
numerator P^. and denominator Q^, 

Ck = ^, (35) 

Qk 

are then defined by the recursions 

Pk^auPk-i+Pk-2 . P-2 = 0, P-i = l, A: = 0,1,2,... (36) 

Qk = akQk-i+Qk-2 , 0-2 = 1, 0_i=0, fc = 0,1,2,... (37) 

Note that these formulas apply equally if fi is rational or irrational; of course 
depending whether fj, is rational or irrational the continued fraction expansion does 
or does not terminate. In the rational case under present discussion we also introduce 
another sequence of nonnegative integers: 

bk^q-H'' [{p-q) Pk-2+qQk-2] , fc = 0,l,2,... (38) 
Given b, let the integer h and the period T{b) be defined by the following formulas: 
bh<b< bh+i, (39) 

T(6) = P,.-2+([^^J +2)P,_1. (40) 

Then the roots of the polynomial can have only one of the following three periods j: 

j = T{b) or j^T{b)~Ph^i or ~] = Ph-i- (41) 

This is the generic case; there are however some cases in which the roots have only the 
two periods T(6) and Pu-i- This happens of course when T{b) = 2Ph-i and whenever 
b takes the following special values: 

b = bh + n{q — bh+i), < n < o/i — 1, n integer. □ (42) 

Remark 4. In case (i) of Proposition 4 the mechanism that yields periods 
longer than unity is the coming into play of the b square-root branch points enclosed 
inside the circle S, which cause a certain number of roots Wjit) to exchange pairwise 
their roles through the time evolution. But this phenomenology only affects some 
roots; others remain unaffected, hence their periods remain unity, and this explains 
the first entry in 1)301) . The precise form of the other entries in this formula, H30() . 
requires of course a more detailed treatment, see [4]; the outcome there depends on 
how many pair exchanges actually do take place, or, equivalently, how many sheets of 
the Riemann surface get actually visited, and this depends in a fairly detailed manner 
on the specific structure of this surface. But note that only two different periods may 
emerge, differing by only one unit. 

In case (ii) of Proposition 4 the second mechanism, associated with the presence 
of the branch point at ^ = 0, comes additionally into play, causing the connection of 
all the b sheets containing the b square-root branch points, both among themselves 
and with the {p — q) sheets containing the branch point at ^ = 0. The corresponding 
(b -\- p — q) roots get permuted among themselves through the time evolution, with 
the period indicated by the second entry in H31() (note that whenever p is quite large 
and b is close to its maximal value q — 1, the resulting period is quite large). The 
remaining (g — 6) sheets, corresponding to the {q — b) branch points lying on the circle 
B outside the circle 5, are isolated, hence the corresponding roots do not take part 
in the quadrille, so that their period remains unity, as indicated by the first entry in 
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The cases (i') and (ii') of Proposition 4 require no additional discussion. 

In case (iii) of Proposition 4 (with < p < q, i. e. < /i < 1) there is no 
branch point at ^ = 0; hence the mechanism is now absent that previously caused 
the connection of all the sheets associated with the b square-root branch points sitting 
on the circle B inside the circle S. This implies, see [4], that each sheet of the 
Riemann surface contains only two branch points, whose projections on the circle 
B are separated by {p — q) other branch points. Note that this entails that two 
branch points which are adjacent on the circle B in the complex ^-plane are instead 
topologically far apart on the Riemann surface of the function w{^), living on sheets 
which are not directly connected. As a consequence the period of the time evolution of 
the root w{t) does not change, as the initial data change causing the circle S to change 
its position and dimension so that the number b of square-root branch points enclosed 
in it increases, until some crucial square-root branch point gets thereby included inside 
the circle S, causing the connection of two separate groups of connected sheets; and 
this mechanism occurs more and more frequently as b increases more and more. This 
explained qualitatively the piecewise constant behavior of the period as b increases, 
characterized by shorter and shorter steps and by bigger and bigger jumps, see (j4()|l - 
(|41|l. The exact treatment of this mechanism, yielding H4()|l - H41|l . is rather complicated, 
as indicated by the role played by the continued fraction expansion of the number — ^ . 
The details are given in [4]. Here we limit ourselves to emphasizing that, for given 
initial data, the generic root wit) can have only 3 possible periods, the third of which 
is just the sum of the first two, see (|41|) .D 

Let us now discuss the case in which /i is an irrational number, recalling that we 
are now considering initial data (|f| — \R\ < |rf,| < |^| + \R\) such that the two circles 
B and S do intersect each other (the results for the other cases have been given in 
Propositions 1-3. ). One must then introduce the ratio v of the length of the arc of 
the circle B that is inside the circle S, to the length of the entire circle B: of course 
< u < 1. Expressing this ratio in terms of the initial data is an exercise in elementary 
plane geometry yielding the formula 

1 



arccos 





2 




2 







(43) 
The 



where the determination of the arccos function must be chosen so that < < 1 
results for the periods are then given by the following 

Proposition 5. (i) If /j, > 1 and the origin ^ = is outside the circle S (i. e 
I77I > 1), the time evolution of the generic root w{t) is still periodic with period T = 
j T and j can take one of the following 3 values: 



i = 1 or j 



1 or 



1 



(44) 



Note that the second and third entry only differ by one unit and moreover that, if 
fj, > 2, the floor functions vanish, hence for all values of /i larger than 2 (and \ri\ > 1) 
the only possible values for j are 1 or 2. 

(ii) If /i > 1 and the origin ^ = is inside the circle S (i. e. \r]\ < 1), the time 
evolution of the generic root w{t) is either periodic with period T or aperiodic. 

(i'). If < and the origin ^ = is outside the circle S (i. e. |?7| > 1), the time 
evolution of the generic root w{t) is still periodic with period T = jT and j can take 
one of the following 3 values: 



1 or j 



1 or j 



(45) 
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Again the second and third entry only differ by one unit and moreover if > 2 the 
floor functions vanish, hence we conclude that for all values of fj, smaller than —2 (and 
\r]\ > 1) the only possible values for j are 1 or 2. 

(ii') If yU < and the origin ^ = is inside the circle S (i. e. \r]\ < 1), the time 
evolution of the generic root w{t) is either periodic with period T or aperiodic. 

(iii). The case < < 1 is again more intriguing, and it requires again the use 
of the continued fraction expansion of jhjii which however now does not terminate. 
We define now in addition the (endless) sequence of real numbers 

iyk = l-i-)^ Pk-2 + Qk-2], fc = 0,l,2,... (46) 

(entailing vq = 0), and we then identify the nonnegative integer h via the inequalities 

0<vh<i^ < Vh+i < 1 • (47) 
Let T{v) be defined by the following expression: 



T{u) = Ph-2 + 



1 - Vh+i 



2 Ph-i- (48) 



Then the motion of the generic root w{t) is again periodic with period T = j T where 
j can take one of the following 3 values: 

j = T{i^) or ~j=T{u)-Ph-i or ~j = Ph-i- □ (49) 

Remark 5. The results for this case with irrational /i can be obtained from those for 
rational /i (see Proposition 4) by taking appropriately the limit in which 

(a) the integers p and q diverge with their ratio fixed (see (|^ ). and 

(b) the number b of square-root branch points inside the circle S, as well as the total 
number q of square-root branch points, also diverge with their ratio fixed (recall 
that all these square-root branch points sit, densely equispaced, on the circle B 
in the complex ^-plane although on different sheets of the Riemann surface of 
the function w (^) - hence this ratio coincides with the quantity v defined, and 
evaluated in terms of the initial data, above, see H43I) '). 

This also suggests obvious extensions to the present case with irrational pi of comments 
contained in the Remark 4, which will not be repeated here. We therefore limit below 
our remarks to aspects of the results reported in Proposition 5 having no immediate 
counterpart in the comments contained in Remark 4. 

In the cases (i) respectively (i') of Proposition 5 the rules giving the period 
of the time evolution of the generic root w(i), see (I44|l respectively (|45|l . are fairly 
straightforward and generally yield rather small periods, unless /i = 1 + e respectively 
fi = —e with e an irrational number positive but extremely small. 

In cases (ii) and (ii') of Proposition 5 the situation is quite interesting because 
the time evolution of the generic root 'w{t) can be either periodic with the basic period 
T or aperiodic. Note that in these cases the circle S intersects the circle B that is 
densely filled with square-root branch points, and moreover the branch point at ^ = 
(which is now of irrational exponent, see ()29|l 'l is inside the circle S. This entails that, 
of the infinity of square-root branch points located on the piece of the circle B that 
is inside S, either none, or all, are active. The first case obtains if the root 'w{t) 
under consideration is initially on a sheet containing a branch point that does not 
fall inside S, hence the time evolution of this root w(t) = w[£,{t)] as the point £,{t) 
travels round and round on the circle S brings it back to its point of departure after 
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a single round; equivalently, in this case the root w{t) = w [^{t)] belongs to a set of 
roots that does not get permuted as the point ^(t) travels round and round on the 
circle S. In the second case the root w{t) under consideration starts from a sheet of 
the Riemann surface that contains a branch point inside B, so that, when ^(t) travels 
round and round on the circle S, an endless sequence of different sheets get accessed 
by w{t) = w [£,{t)]; equivalently, such a root w{t) = w [£,{t)] belongs to a set (including 
an infinity of roots) that does get permuted as the point ^ (t) travels round and round 
on the circle S, with both mechanisms - the pairwise exchange of some roots, and 
the cyclic permutation of an infinite number of roots - playing a role at each round. 
The identification of which sheets get thereby accessed, and in which order - namely 
the specific shape of the trajectory when looked at, as it were stroboscopically, at 
the discrete sequence of instants Tk = kT, k = 1,2,3,...- is discussed in [4]. The 
extent to which this regime yields irregular motions is discussed further below, also 
to illuminate the distinction in these regimes between the time evolution entailed by 
our model with a given irrational value of fi, and that of the analogous models with 
rational values of fi providing more and more accurate approximations of the given 
irrational value of /i. 

In case (iii) the time evolution is still isochronous, inasmuch as the results reported 
above entail that, for any given initial data (excluding, of course, the special ones 
leading to a collision; which are special in the same sense as a rational number is 
special in the context of real numbers), the motion of every root 'Wj{t) is periodic 
with one of the 3 periods entailed by (and note that the value of the integer 
j provided by the third of these 3 formulas is just the sum of the 2 values for j 
provided by the first 2 of these 3 formulas). It is indeed clear that the initial data 
yielding such an outcome are included in an open set of such data, having of course 
full dimensionality in the space of initial data, all yielding the same outcome: since 
the periods do not change, see (|49|) . if the change of the initial data, hence the change 
in the ratio is sufficiently tiny. However the measures of these sets of data yielding 
the same outcome gets progressively smaller as the predicted periods get larger, and 
moreover the corresponding predictions involve more and more terms in the (never 
ending) continued fraction expansion of the irrational number jz^, see (|34f) . displaying 
thereby, as u increases towards unity, a progressively more sensitive dependence of the 
periodicity of our system on the initial data and moreover on the parameters (the 
coupling constants, that determine the value of fi, see CU) of our physical model 

Example. Let us display here a specific example with the following (conveniently 
chosen) irrational value of fi in the interval < /i < 1 (hence corresponding to case 
(iii) of Proposition 5): 

2 1 1 1 + V5 

3 + v5 1 + ip 1 - fi 2 
where of course ip is the golden ratio. This assignment entails that all the coefficients 
ak, see (|34|) . are in this case unity, = 1, hence the quantities Pk, see (|37|) . coincide 
with the Fibonacci numbers F^, 

Pk = Fk , F_2 = , ^^-1 - 1 , -Fo = 1 , (51) 

Fk+2 = Fk + Fk+1 , (52) 
and moreover one easily finds that 

i^k^l-ip-''- (53) 
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The corresponding formula for the periods obtains inserting these values in H47|) and 
H49|l . From it, with some labor, one can obtain the following controlled estimate for 
the possible values of j: 

„{2-i.) V5 + 2-(75 + l) i.i2-u) ^^^^ 




V5(l-i^)-'- (75-1) (1- I/) 

These inequalities are valid for all values of i/ (in the interval < < 1); they clearly 

entail that the integer j diverges proportionally to (1 — i^)^^ as v 1, and that for 
< v < P with 

-- 0.760067 . . . (55) 

the only possible values of j are 2 and S.D 

This concludes our presentation of the results, and of some related observations, 
detailing the periodicity (if any) of the time evolution of a generic root wit) of (|16|l 
with H15|l . The identification of analogous, but of course more definite, results for the 
physical root w{t), and the consequential information on the periodicity (if any) of 
the solution of the physical problem - as well as some additional information on 
the corresponding trajectories of the coordinates Zn{t) - are provided in [4]. 

Last but not least let us elaborate on the character of the aperiodic time evolution 
indicated under item (ii) of Proposition 5 , including the extent it is irregular and it 
depends sensitively on its initial data. It is illuminating to relate this question with 
the finding reported under item (ii) of Proposition 4, also in order to provide a better 
understanding of the relationship among the aperiodic time evolution that can emerge 
when fjL is irrational (see item (ii) of Proposition 5) and the corresponding behavior 

- say, with the same initial data - for a sequence of models with rational values of 
/z (see (Unj) that provide better and better approximations to that irrational value 
of /i; keeping in mind the qualitative difference among the aperiodic time evolution 
emerging when /i is irrational, and the periodic - indeed, even isochronous - time 
evolutions prevailing whenever (x is rational, albeit with the qualifications indicated 
under item (ii) of Remark 4. Note that we are now discussing the case /z > 1 (an 
analogous discussion in the < case can be forsaken), with initial data such that 
the two circles B and S in the complex ^-plane do intersect and moreover the origin 
^ = falls inside the circle S (i. e. |rf,| < |f | + \R\ and \ri\ < 1). 

Let us then consider a given irrational value of > 1 and let the rational number 

- (with p > q) provide a very good approximation to /Lt, which of course entails that 
the positive integers p and q are both very large. Consider then the difference 

Aj = A6 (56) 

(see the second entry in ()31|l) between the two positive integers that characterize the 
two periods of the two time evolutions of w{t) corresponding to two sets of initial 
data that differ very little. Here clearly the quantity A 6 is the difference between the 
number of branch points that are enclosed inside the circle S for these two different 
sets of initial data. Since the number q of branch points on B is very large, it stands 
to reason that 

Ab^lOiqS)\ , (57) 

where the quite small (positive) number S provides a (dimensionless) measure of the 
difference between the two sets of initial data (see for instance (j^ . which clearly 
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becomes approximately applicable when q is very large). The floor symbol [J has 
been introduced in the right-hand side of this formula to account for the integer 
character of the numbers b hence of their difference A 6, while the order of magnitude 
symbol O ( ) indicates that the difference A 6 is proportional (in fact equal, given the 
latitude left by our definition of the quantity S) to the quantity q 5 up to corrections 
which become negligibly small when S is very small and q is very large, but irrespective 
of the value of the quantity q 5 itself which, as the product of the large number q by 
the small number S, is required to be neither small nor large. 
The relation 

Aj=lO{qd)\ (58) 

implied by this argument indicates that, for any given q, one can always choose (finitely 
different) initial data which differ by such a tiny amount that the corresponding periods 
are identical, confirming our previous statement about the isochronous character of 
our model whenever the parameter /i is rational. But conversely this finding also 
implies that, for any set of initial data in the sector under present consideration (i. e. 
that characterized by the inequalities |?7| < 1 and |f | + > |rf,| , and by an additional 
specification to identify the physical root w{t), see [4]), if our physical model (Q) is 
characterized by an irrational value of /i, see (|17|l . and one replaces this value by a 
more and more accurate rational approximation of it, see (|20|l ~ as it would for instance 
be inevitable in any numerical simulation - corresponding to larger and larger values 
of p and q, then one shall have to choose the two different assignments of initial data 
closer and closer to avoid a drastic change of period - and for these very close sets 
of data the motion is indeed periodic with a period (which we are able to predict, 
see H31|) and [4], but) which becomes larger and larger the better one approximates 
the actual, irrational value of fi. Moreover in any numerical simulation the accuracy 
of the computation, in order to get the correct period, shall also have to increase 
more and more (with no limit), because of the occurrence of closer and closer near 
misses through the time evolution (associated with the coming into play of active 
branch points sitting on the circle B closer and closer to the points of intersection 
with the circle S). And finally, if one insists in treating the problem with a truly 
irrational pL, then, no matter how close the initial data are, the change in the periods 
becomes infinite because the difference A6 in the number of active square-roots branch 
points on the circle B included inside the circle S is infinite (see (ISSf) ). signifying that 
the motion is aperiodic, and that its evolution is indeed characterized by an infinite 
number of near misses, making it truly irregular. 

This phenomenology, together with that of the near misses as described above, 
illustrates rather clearly the irregular character of the motions of our physical 
model when the coupling constants have appropriate values (such as to produce an 
irrational value of outside the interval < /i < 1) and the initial data are in 
the sector identified above. Note that the Lyapunov coefhcients associated with the 
corresponding trajectories vanish, because these coefficients - as usually defined - 
compare the difference (after an infinitely long time) of two trajectories that, to begin 
with, differ infinitesimally; whereas our mechanism causing the irregular character of 
the motion requires, to come into play, an arbitrarily small but finite difference among 
the initial data. The difference between these two notions corresponds to the fact that 
inside the interval between two different real numbers - however close they may be - 
there always is an infinity of rational numbers; while this is not the case between two 
real numbers that differ only infinitesimally! This observation suggests that, in an 
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applicative context, the mechanism causing a sensitive dependence on the initial data 
manifested by our model may be phenomenologically relevant even when no Lyapunov 
coefficient, defined in the usual manner, is positive. As already observed previously [2], 
this mechanism is in some sense analogous to that yielding aperiodic trajectories in a 
triangular billiard with irrational angles; although in that case - in contrast to ours - 
this outcome is mainly attributable to the essentially singular character of the corners, 
and moreover no truly irregular motions emerge. 



3. Outlook 

In this paper we have introduced and discussed a 3-body problem in the plane 
suitable to illustrate a mechanism of transition from regular to irregular motions. 
This model is the simplest one we managed to manufacture for this purpose. Its 
simplicity permitted us to discuss in considerable detail the mathematical structure 
underlining this phenomenology: this machinery cannot however be too simple since it 
must capture (at least some of) the subtleties associated with the onset of an irregular 
behavior. Therefore in this short paper we were only able to report our main findings 
without detailing their proofs, and we also omitted several other relevant aspects of our 
treatment (including a fuller discussion of previous work by others in related areas): 
this material shall be presented in a separate, much longer, paper [4], and probably 
as well via an electronic version of our findings so as to supplement their presentation 
with various animations illustrating these results and their derivation. 

Our main motivation to undertake this research project is the hunch that this 
mechanism of transition have a fairly general validity and be relevant in interesting 
applicative contexts. Hence we plan to pursue this study by focussing on other 
cases where this mechanism is known to play a key role, including examples (see, 
for instance, [2] and [7]) featuring a pattern of branch points covering densely an area 
of the complex plane of the independent variable rather than being confined just to 
reside densely on a line as is the case in the model treated herein; and eventually to 
extend the application of this approach to problems of direct applicative interest. 

In this connection the following final observation is perhaps relevant. In this 
paper as well as in others [1, 2, 7] the main focus has been on models featuring 
a transition from an isochronous to an irregular regime, and in this context much 
emphasis was put on the "trick" ^ and in particular on the relationship it entails 
between the periodicity of the ( "physical" ) dependent variables z„ (t) as functions of the 
real independent variable t ("time") and the analyticity of other, related ("auxiliary") 
dependent variables C„ (r) as functions of a complex independent variable r. But 
our findings can also be interpreted directly in terms of the analytic properties of the 
physical dependent variables z„ {t) as functions of the independent variable t considered 
itself as a complex variable. Then the time evolution, which corresponded to a uniform 
travel round and round on the circle C in the complex r-plane or equivalently on the 
circle S in the complex ^-plane, is represented as a uniform travel to the right along 
the real axis in the complex i-planc, while, via the relations (see lO and (|15(l ) 



the pattern of branch points in the complex r-plane or equivalently in the complex ^- 
plane gets mapped into a somewhat analogous pattern in the complex i-plane, repeated 




(59) 
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periodically in the real direction with period T, see (j^J. In particular - to mention 
the main features relevant to our treatment, see above - the circle B on which the 
square-root branch points in the complex ^-plane sit, gets mapped in the complex t- 
plane into a curve B on which sit the square-root branch points in the complex t-plane; 
note that this curve B (in contrast to the circle B) does now depend on the initial 
data. This curve is of course repeated periodically; it is closed and contained in each 
vertical slab of width T (see figure Et) if the point ^ is outside B , otherwise it is open, 
starting in one slab and ending in the adjoining slab at a point shifted by the amount 
T; and it does not or does cross (of course twice in each period) the real axis in the 
complex <-plane depending whether, in the complex ^-plane, the two circles B and S 
do not or do intersect each other (see figure Likewise, depending whether it is 
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Figure 2. Evolutionary paths and branch points of the sohitions in the complex 
5-plane and t-plane for /i = 15/37. 



inside or outside the circle S, the point ^ = - which, as entailed by our analysis, 
is a highly relevant branch point in the complex ^-plane (unless < < 1) ~ gets 
mapped into an analogous branch point located in each vertical slab above or below 
the real axis in the complex t-plane; while the other branch point, at ^ = oo in the 
complex ^-plane, gets mapped into an analogous branch point located at infinity in 
the lower half of the complex t-plane. Clearly the physical mechanism of near misses, 
which is the main cause of the eventual irregularity of the motion, becomes relevant 
only for initial data such that the curve B crosses the real axis, thereby causing (if fi 
is irrational) an infinity of square-root branch points of the functions Zn{t) to occur 
arbitrarily close to the real axis in the complex i-plane - branch points which are 
however active (namely, they actually cause a near miss in the physical evolution) 
in only some (yet still an infinity) of the infinite number of vertical slabs in which 
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the complex t-plane gets now naturally partitioned. The near miss implies that the 
two particles involved in it slide past each other from one side or the other depending 
whether the corresponding branch point is just above or just below the real axis in 
the complex f-plane. The sensitive dependence on the initial data is due to the fact 
that any tiny change of them causes some active branch point in the complex t-planc 
which is very close to real axis to cross over from one side of it to the other, thereby 
drastically changing the outcome of the corresponding near miss. 

This terse discussion shows clearly that the explanation of the irregular behavior 
of a dynamical system in terms of travel on a Riemann surface is by no means restricted 
to isochronous systems. We found it convenient to illustrate in detail this paradigm 
by focussing in this paper on a simple isochronous model and by using firstly r and 
then ^ as independent complex variables - but, as outlined just above, our analysis can 
also be done - albeit less neatly - by using directly the independent complex variable 
t; and the occurrence of a kind of periodic partition of the complex t-plane into an 
infinite sequence of vertical slabs - characteristic of our isochronous model - does 
not play an essential role to explain the irregular character of the motion when such a 
phenomenology does indeed emerge. The essential point is the possibility to reinterpret 
the time evolution as travel on a Riemann surface, the structure of which is sufficiently 
complicated to cause an irregular motion featuring a sensitive dependence on its initial 
data. The essential feature causing such an outcome is the presence of an infinity of 
branch points arbitrarily close to the real axis in the complex t-plane, the positions 
of which, as well as the identification of which of them are active, depends on the 
initial data nontrivially. The model treated in this paper shows that such a structure 
can be complicated enough to cause an irregular motion, yet amenable to a simple 
mathematical description yielding a rather detailed understanding of this motion; this 
suggests the efficacy also in more general contexts of this paradigm to understand 
(certain) irregular motions featuring a sensitive dependence on their initial data and 
possibly even to predict their behavior to the extent such a paradoxical achievement 
(predicting the unpredictable!) can at all be feasible. 
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